Interacting Dark Energy — constraints and degeneracies 
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In standard cosmologies, dark energy interacts only gravitationally with dark matter. There 
could be a non-gravitational interaction in the dark sector, leading to changes in the effective DE 
equation of state, in the redshift dependence of the DM density and in structure formation. We 
use CMB, BAO and SNIa data to constrain a model where the energy transfer in the dark sector is 
proportional to the DE density. There are two subclasses, defined by the vanishing of momentum 
transfer either in the DM or the DE frame. We conduct a Markov-Chain Monte-Carlo analysis to 
obtain best-fit parameters. The background evolution allows large interaction strengths, and the 
constraints from CMB anisotropies are weak. The growth of DM density perturbations is much 
more sensitive to the interaction, and can deviate strongly from the standard case. However, the 
deviations are degenerate with galaxy bias and thus more difficult to constrain. Interestingly, the 
ISW signature is suppressed since the non-standard background evolution can compensate for high 
growth rates. We also discuss the partial degeneracy between interacting DE and modified gravity, 
and how this can be broken. 
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I. INTRODUCTION 

The late-time acceleration of the expansion of the 
Universe demands explanation and observational veri- 
fication. Currently observational tests of the standard 
ACDM cosmological model are not precise enough to ad- 
equately rule out the wide variety of alternative dark en- 
ergy (DE) models that have been proposed to explain 
the data. Instead it is necessary to obtain constraints 
on the free parameters of such models and find ways to 
distinguish between them using observations. One way 
to test ACDM is to describe DE as an effective fluid and 
promote its equation of state w to a free parameter. This 
is known as u>CDM and allows deviations from the stan- 
dard value of w = — 1; w is still only measured to about 



acy 



some datasets 

Alternatively, one may go beyond wCDM by introduc- 
ing a new parameter to quantify interactions within the 
dark sector. It is natural to expect some new physics in 
the dark sector given the richness of interactions between 
species in the standard model of particle physics [2]. In- 
deed taking dark sector interactions to be zero would be 
an additional assumption of the model. Models of dark 
sector interactions (see e.g. 0-[7^) have little concrete 
guidance from particle physics, but by studying possi- 
ble interactions and confronting them with observations 
we can shed light on questions such as which models may 
lead to unphysical behaviour and which are in best agree- 
ment with observations. 

Further motivation for interacting dark energy (IDE) 
models includes: (1) they may alleviate the coincidence 
problem of explaining why the domination of DE roughly 
coincides with the formation of large-scale structure; 
(2) IDE affects structure formation and therefore pro- 
vides a new way to modify the predictions of the standard 
non-interacting model. 



Here we investigate a p articular version of the inter- 
action model used in |77H85|. The general form of this 
interaction in the background is Tp, where T is constant. 
This has also been used to describe particle decays in 
other contexts (86I-I88}. Interactions of the form aHp 
(with a constant) (see e.g. [89] and references therein) 
appear similar, but they mean that the interaction at 
any event is affected by the global expansion rate H . as 
opposed to the locally determined interaction Tp. 

We describe the IDE model in the background uni- 
verse in Sec. HTland in the perturbed universe in Sec. IIIII 
In Sec. IIVI we investigate the effects of the interaction 
on the CMB and matter power spectra. We find the 
best- fit models using CMB, baryon acoustic oscillation 
(BAO) and supernova (SNIa) data, and we discuss the 
behaviour of some typical models and the implications 
for the growth of large-scale structure in the Universe. 
Our conclusions are in Sec. [V] 



II. IDE IN THE BACKGROUND 

An IDE model is characterized in the background by 
the energy transfer rate Q x = — Q c : 



p' c = -3Hp c 
p' x = -3H(1 



aQ c , 
w)p x 



aQ x , 



(1) 
(2) 



where primes denote derivatives with respect to confor- 
mal time t,T-L = dlna/dr and w = p x /p x . We can define 
an effective DE equation of state parameter to be that of 
a non-interacting DE with the same p x (a), i.e. 



aQ x 
"H-Lpx 



(3) 



We can see from Eq ([3]) that w e {£ can be dynamical even 
if w is a constant. Interestingly w e & can be less than — 1, 
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or cross —1 during its evolution if Q x > 0, even though 
w itself is always greater than — 1. 

A simple model for Q x is a linear function in the dark 
sector energy densities. IDE with Q x oc p c has been 
studied in the greatest detail [77M84l |. However, for con- 
stant w the model suffers from an instability [79]. This 
instability arises because the model of DE as a fluid with 
constant w is non-adiabatic. The instability can be cured 
by allowing w to vary in time [83| . 

Here we study the version with Q x oc p x , 



Q, 



-Qc = Fp x , 



(4) 



where T is a constant transfer rate. The strength of the 
interaction is measured by |r|/iJ - T > corresponds to 
energy transfer from DM— >-DE. This appears somewhat 
unnatural, since the energy transfer is proportional to p x . 
For r < 0, the interaction can be seen in the background 
as a decay of DE into DM, which is a more natural model. 
The solution of © is [13 



Px — Pxoa 



-3(l+u 



'exp [r(t-*o)], 



(5) 



which shows that T > leads to exponential growth of 
DE. By ([lj) , it follows that p c eventually becomes nega- 
tive. The model breaks down if this happens before the 
current time, which is possible for large T/Hq- Observa- 
tional constraints require T / Hq < 1, so that typically the 
DM density only becomes negative in the future. In this 
case, we can treat the model as a viable approximation, 
for the past history of the Universe, to some more com- 
plicated interaction that avoids the blow-up of DE in the 
future. The DE-+DM decay model, with T < 0, does not 
have this problem: both energy densities remain positive 
at all times when evolving forward from physical initial 
conditions [80]. Furthermore, the T < case includes the 
possibility of beginning with no DM present and having 
it created entirely from the decay of DE. 

We use a phenomenological fluid model for DE, in 
which we treat w and the soundspeed c s as arbitrary 
parameters. This is a commonly used model for non- 
interacting DE, where the model is known as uiCDM. We 
impose the condition w > —1 to avoid 'phantom' insta- 
bilities that can arise in scalar field models of DE [90l.[9l|. 
The limiting case w = — 1 is admitted by the background 
equations, but the perturbation equations have singular- 
ities (see below). Therefore we assume 



w > —1, w = const. 



(0) 



For completeness, we consider also the w < — 1 case in 
Appendix [A"l In the background, the T < case appears 
to be better motivated. However, the analysis of pertur- 
bations (see below) shows that these models suffer from 
an instability when w > — 1. The T > models avoid 
this instability. 



III. IDE IN THE PERTURBED UNIVERSE 

The critical difference between the background and 
perturbed IDE is that there is nonzero momentum trans- 
fer in the perturbed universe. As emphasized in [79], a 
model for energy and momentum transfer does not follow 
from the background model - and a covariant and gauge- 
invariant approach is essential to construct a physically 
consistent model for energy-momentum transfer. 

General IDE 

We give a brief summary of the general discussion 
in [791 ]. The perturbed Friedmann metric in a general 
gauge is 

ds 2 = a 2 { - (l + 2(f)) dr 2 + 2d l Bdrdx l 

+ [(1 - 2V>)% + 2d i d j E]dx i dx j }, (7) 

Each fluid A satisfies an energy-momentum balance equa- 
tion, 

v„r2" = Q% X A Q A = 0, (8) 

where the second condition expresses conservation of the 
total energy-momentum tensor. For dark sector interac- 
tions, the energy-momentum transfer four- vectors satisfy 

= -Qc- (9) 

We split Q^ A relative to the total four- velocity w M , so that 
Qa - Qau»+F a \ Q a = Qa+SQa, u^F% = 0, (10) 



where Qa is the energy density transfer rate relative to 

A 



and F A is the momentum density transfer rate relative 



to u M . To first order 

F A i = a- 1 (0,d i f A ), 



(11) 



where ] a is the (gauge-invariant) momentum transfer po- 
tential. 

We choose each u A and the total it^ as the unique 
four-velocity with zero momentum density, i.e. 

n u u A = -PAU A , TW = -pu», (12) 

PA = PA + SpA, p = ^APA = P + Sp- (13) 

Then we have 



a' 1 (1 - <j), d l v A ) , = a' 1 (l - tf>, d l v) , (14) 



(T, A PA + T, a pa)v = Y> a (pa +pa)va, 



(15) 



where va,v are the peculiar velocity potentials. Equa- 
tions ((10]) and flU) imply that 



1 



6Qj 



Qf = ad, Q A (v + B)+f A 



(16) 
(17) 
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The evolution equations for S A and the velocity per- 
turbation 



= -k 2 (v A + B), 



(18) 



are |7E 



S' A + 3U(c 2 sA - w A )8 A + (1 + w A )6 A 
+ 9n 2 (l + w A )(c 2 sA ~c 2 aA )^ 
- 3(1 + w A )ip' + (1 + w A )k 2 (B - E') 



clQa 

PA 



5 A + m{c 2 s 



+ —SQ A , (19) 

PA 



0' A + H(l-3c 2 sA )9 A 



-^^k 2 5 A -k 2 , 
(1 + w A ) 



(1 + w A )p J 



-{Qa[o-( 



i + c: 



a)Va\ 



2 /a}- (20) 



Here c sA is the physical soundspeed, defined by c 2 A = 
(5p A / 5 ^PA)restframo, and c aA is the adiabatic soundspeed, 
defined by c 2 A = p' A /p' A - For the adiabatic DM fluid, 
c 2 c = c 2 c = w c = 0. By contrast, the DE fluid is non- 
adiabatic: c 2 x = w < and so c ax cannot be the physical 
soundspeed. The physical soundspeed for the fluid DE 
model is a phenomenological parameter. It must be real 
and non-negative to avoid unphysical instabilities. We 
choose c sx — 1, which is the sounspeed for quintessence (a 
self-consistent model of DE). Our analysis is insensitive 
to the value of c sx , as long as c sx is close to one, so that 
DE does not cluster significantly on sub-Hubble scales. 
(See [79( for more details.) 



DM-baryon bias from IDE 

In IDE models, the DE exerts a drag on DM but not 
on baryons. This leads to a linear DM-baryon bias in the 
late-time density perturbations, and in general also to a 
velocity difference [82]. For baryons after decoupling 

8' b +0 b -3ip' + k 2 {B-E') = o, e' b +ne b -k 2 (j) = o. (21) 

Thus for non-interacting DE models, 



- 9b — (&c - 9b)i — 
a 



We can choose (9 b — c )i = 0, so that 

e c - e b = o, 6 c -s b = (6 C - <5 b )i 



(22) 



(23) 



Thus in standard DE models, there is no DM-baryon 
velocity difference, and any linear density perturbation 
difference is determined by initial conditions. 

For IDE models, the interaction induces a non- 
constant difference between 8 C and 8 b - which is degen- 
erate with the standard galaxy bias. The Euler equation 



for DM is 0, with c 2 s 



Wr 



0. This differs from the 



standard Euler equation unless k 2 f c = Q(9 — 8 C ), which 
follows only for Q£ — Q c u^, regardless of the form of 
Qc [13] • In those models that modify the Euler equation 
for DM, there will also be a velocity bias. Equations (QI 
and ([H imply 



(8 c -5 b y + ^(S c -6 b ) + (< 

Pc 

"' \SQ C + Q C (4>- §b) 



Pc 



a 

Pc 



k 2 f c + Qc(9 - 



(24) 



(25) 



Thus there will be a velocity bias, unless = Q 



Our IDE models 

The preceding equations are completely general. A 
choice must now be made for the energy-momentum 
transfer in the dark sector. Firstly, the nature of the 
background energy transfer suggests that we take 



Qx = ^Px = r j o x (l + S x ) 



-Qc 



(26) 



where 8 A = 8p A / p A . Thus we are treating T as a univer- 
sal constant. For the momentum transfer, the simplest 
physical choice is that there is no momentum transfer in 
the rest frame of either DM or DE (zH This leads 
to two types of model, with energy-momentum transfer 
four-vectors parallel to either the DM or the DE four- 
velocity: 



Qx 



QxK = type: Q||u c 

Qx< = -QV type: Q|| Uj 



Thus 



= aTp x [l + 8 x + 0, di(v A + B)] . 



(27) 
(28) 



(29) 



where A — c, x for type Q\\u c , Q\\u x . By (JT7J) , the 
momentum transfer relative to the background frame is 

fx = Tp x (v c - v) = —f c for Q\\u c , (30) 
fx = Tp x (v x - v) = -f c for Q\\u x . (31) 

For both the Q||u c and QHw^ models, the density per- 
turbation (continuity) equation (|19p reduces to 

8' c + 6 C - 3V>' + k 2 (B - E') = aT^(8 c - 8 X - <j>), (32) 

Pc 

Q 

8' x + m(l - w)8 x + (1 + w)6 x + 9H 2 (l - w 2 )^ 



- 3(1 + w)ip' + (1 + w)k 2 (B - E 1 ) 
= aT U + 3H(l-w)?| 



(33) 
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The velocity perturbation (Euler) equations are however 
different. For the Q||w c model, (f2"0")) gives 



9' c + ne c - k 2 <f> = o, 



(1+w) (1 + w) 

For the Q\\u x model: 



(34) 

- 20*). (35) 



9' C + W C -k 2 4> = aT^(9 c -e x ), (36) 



6' x - 2W X 



k 2 5 T 



(1 + w) 



Pc 

aT9 x 
'(1 + w)' 



(37) 



It follows that the Euler equation for DM in the 
Q||u c model has the standard form, whereas it is modified 
in the Q||u x model. 



Instability 

There is an obvious instability in the Euler equations 
for DE, (|3"5|) and (|37p . as w — > — 1. Thus we must exclude 
the value w = — 1. This instability is different from that 
for a dynamical DE model with w crossing —1, in which 
case the DE perturbation is well-defined, but at least one 
more degree of freedom is required, usually leading to its 
interpretation as a sign of modified gravity. Here though, 
the DE is not dynamical and the DE perturbation is ill- 
defined at w = — 1. 

These equations also reveal an instability for w ^ — 1 
in certain regions of parameter space. The underlying 
cause of this instability is the choice of c 2 x — 1, which 
means that the DE fluid is non-adiabatic, as discussed 
above. It is qualitatively similar to the instability first 
discovered for constant w IDE in [7^|. (See also (34. l92l- 
|96| for the case of models with T replaced by aH) . This is 
a DE velocity instability, which then drives an instability 
in the DE and DM density perturbations. 

On large scales, we can drop the 8 X and <f> terms in the 
DE Euler equations (P5"j) and (|37[) . In (p?5j) we can also 
set C = by (fM| . Then we can integrate to find that 



a(r=o) 



exp 



1 + w 



(38) 



where 6 X T °' is the DE velocity in the non-interacting 
case, and a = 2,1 for Q\\u c , Q\\u x . It follows that 

r 



(1 + w) 



> 



instability. 



(39) 



Note that although one can choose a reference frame 
where 9 X = 0, the instability is still present in the ve- 
locity difference, which is gauge invariant. Given our 
assumption ©, the stable models must have positive T, 
i.e. 



for both Q\\u c and Q\\u x . This defines for us the phys- 
ically acceptable models. In Appendix [A] we allow for 
any sign of T and 1 + w. In order for the instability to 
affect significantly the perturbation evolution by today, 
the time scale of growth of 9 X in (|38|) should be shorter 
than the Hubble time, i.e., the models with 



> 



H {l + w) 



1 for Q||u c , 

2 for Q\\u x . 



(41) 



w > — 1 and T > 



no instability, 



(40) 



may not be viable. The results from our full parameter 
scan confirm this (see the excluded wedges near to w = 
-1 in Fig. U). 



IV. ANALYSIS 

The evolution of TwCDM models was computed nu- 
merically using a modified version of the CAMB Boltz- 
mann code [97| . including implementation of the initial 
conditions derived in Appendix [B] The code was adapted 
(a) to allow for the non-standard background evolution 
caused by the interactions; (b) to evolve the DM veloc- 
ity perturbation (ordinarily set to zero); (c) to suppress 
perturbations when |1 + w\ < 0.01 due to the blow up of 
terms in ([33)1 and (pT7f as w — > — 1. It is useful to include 
the w = —1 limit for comparison with ACDM. 

Insight into the physical implications of the interaction 
can be gained by running the modified CAMB code with 
fixed input parameters, varying only the interaction rate 
r. Figure [1] shows the CMB power spectrum for three 
values of T with all other cosmological parameters set to 
typical values (see Table Q] for details). 

Positive T describes a transfer of energy from DM 
to DE, so with fixed il c today, the DM energy density 
would have been correspondingly greater in the past than 
without interactions. Hence the amplitude of the CMB 
power spectrum is decreased and the position of the peaks 
shifted, since a larger proportion of DM at early times 
implies a smaller amount of baryonic matter and there- 
fore a more significant effect from photon driving before 
decoupling. The present-day matter power spectrum for 
these choices of T shows that a relative increase in the 
past DM density naturally leads to more structure for- 
mation and an increase in the amplitude of the matter 
power spectrum. 

The modified CAMB code was integrated into the Cos- 
moMC HI Markov-Chain Monte-Carlo (MCMC) code 
in order to explore the parameter space (see Appendix IA1 
for more details). The data used in the MCMC analy- 
sis were: CMBQVMAP7 Q), BAO HST 03, and 
SNIa (SDSS [10l|) dat a, as well as a prior on fif, from big- 
bang nucleosynthesis |102| . Figure [5] shows the 68% and 
95% likelihood contours in the w — T plane for the two 
different momentum transfer models Q||u c and , 
where all other parameters have been marginalized over. 

The likelihood regions are very similar for the two 
models since they differ only in their perturbations 



5 



Model 


^ A 


Ay 2 


T/H 


w 


H 


Q. b h 2 


Q c h 2 


n s 


A s 




ACDM best-fit 









-1 


69.8 


0.0223 


0.113 


0.960 


2.16xl0 -9 


0.0844 


ACDM69 




0.77 '4 




-1 


69.0 


0.0221 


0.114 


0.958 


2.18xl0~ 9 


0.0855 


ACDM70 




-0.0200 




-1 


70.0 


0.0224 


0.112 


0.962 


2.16xl0 -9 


0.0844 


ujCDM best-fit 




-0.220 




-1.03 


70.7 


0.0222 


0.113 


0.960 


2.18xl0~ 9 


0.0883 


TiuCDM A 


Q \\u c 







-0.98 


70.0 


0.0226 


0.112 


0.960 


2.10xl0~ 9 


0.0900 


TtoCDM b 


Uc 




0.2 


—0.98 


70.0 


0.0226 


0.112 


0.960 


2.10xl0~ 9 


0.0900 


TwCDM C 


Q\\u c 




0.4 


-0.98 


70.0 


0.0226 


0.112 


0.960 


2.10xl0~ 9 


0.0900 


TwCDM la 


Uc 


-0.00830 


0.4 


—0.95 


70.9 


0.0222 


0.0702 


0.961 


2.16xl0~ 9 


0.0816 


ItoCDM lb 


Q\\u c 


0.702 


0.7 


-0.85 


70.0 


0.0223 


0.0311 


0.963 


2.15xl0" 9 


0.0832 


TwCDM 2a 


Q \\u x 


-0.236 


0.4 


-0.95 


71.0 


0.0224 


0.0701 


0.966 


2.19xl0 -9 


0.0870 


TiuCDM 2b 


Q\\u x 


-0.0420 


0.7 


-0.85 


70.2 


0.0224 


0.0305 


0.966 


2.15xl0" 9 


0.0819 


r > 0, w > -1 best-fit 


Q Me 


-0.0522 


0.366 


-0.964 


71.0 


0.0224 


0.0748 


0.963 


2.18xl0~ 9 


0.0849 


r > 0, w > -1 best-fit 


Q\\u x 


-0.322 


0.798 


-0.851 


70.4 


0.0224 


0.0194 


0.965 


2.18xl0 -9 


0.0870 



TABLE I: Cosmological parameters for IDE models (see Appendix [A] for more general constraints). 
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FIG. 1: CMB and total matter power spectra from the modified CAMB code for 3 Q||u c models with different values of Y but 
identical values of their remaining parameters (see FwCDM A,B,C in Table[T]). 




FIG. 2: Smoothed 68% and 95% contours of the marginalised probability distribution for IDE model with Q||u c (left) and 
(right) in the range of stability, w > —1 and V > 0. Crosses identify models chosen to be analyzed in more detail (see 
Table I). 
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FIG. 3: CMB and total matter power spectra from the modified CAMB code for the WMAP7 iuCDM best-fit values and 
the FujCDM la,lb,2a,2b models chosen from the 95% confidence range for further analysis (see Table [I]). The best-fit values 
of standard cosmological parameters were found using CosmoMC. Models la,lb have F = 0.4ffo and Q\\u c while 2a, 2b have 
T = Q.THq and Q\\u x . 




FIG. 4: Left: Normalized growth rates for ACDM and the same best-fit models as in Fig. [3] Right: The same models but 
showing a normalized combination of a 2 p c S c which is important for the ISW effect. 
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FIG. 5: Deviations from ACDM of the effective Hubble parameter (left) and effective Newton constant for S c (right), for the 
same best-fit models as in Fig. 



7 



2 4 6 8 10 



-0.85 



-0.90 



-0.95 



eff 



■1.05 



1.10 



TwCDM lb,2b 



TwCDM la,2a 



_ 3 
p a 



ACDM 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 



/ / 
// 

/ / 
/ / 



ACDM 
TwCDM la,2a 
TwCDM lb,2b 



FIG. 6: Comparison with ACDM of the effective DE equation of state (left) and the DM density (right), for the same best-fit 
models as in Fig. [3] 



and the observations predominantly constrain the back- 
ground evolution. The best- fit values for the Q||m c and 
Q\\u x models are different due to the Integrated Sachs- 
Wolfe (ISW) effect on the CMB, as shown in Table HI 
For the Q\\u x model the best-fit value is a genuine global 
maximum. For the Q||m c model however the mean like- 
lihood function of w and T is essentially one-tailed, with 
the true global maximum lying outside of the region we 
consider physical (see Appendix Indeed the \ 2 of this 
point is close to that of ACDM (see Table I). This is be- 
cause <2||ii c models in this region can closely mimic the 
ISW signature of ACDM. 

The CMB data is best fit by a particular ISW sig- 
nal, and the two momentum transfer models differ some- 
what in their structure formation histories. DM in the 
<5||u x model receives a change in momentum from the DE 
perturbations, as expressed by its modified Euler equa- 
tion (|34f . leading to more structure growth relative to 
the Q||it c model. This means that DE can be weaker for 
the Q\\u x model in order to give the same amount of ISW 
signal as the Q\\u c model. 

In order to assess the relative merits of these models 
we have included the change in x 2 from a ACDM base- 
line. To help put this quantity into context we have also 
included two best-fit ACDM models with Hq fixed at 69 
and 70 km/s/Mpc. The mean likelihoods of the sam- 
ples vary little in the direction of the degeneracy in the 
w — T plane. For example, the difference in A% 2 between 
the Q\\u x best-fit and the ACDM best-fit is less than 
the difference between the two fixed- Hq ACDM models 
(ACDM69 and ACDM70 in Table 0J. 

In Fig. 6, we show the effective DE equation of state 
and the a 3 -scaled energy density for DM for a selection 
of models in comparison with ACDM. Interestingly, we 
find that the w e g for the best-fit TwCDM models with 
w = —0.85 and —0.95 crosses —1 during its evolution, 



showing a quintom-like behaviour [lQ3j- 

We have focused on the stable T > models with 
w > — f . These models do have a problem of negative 
DM energy densities in the future, but we assume that 
this can be cured by a more realistic model to which our 
model is a good approximation when p c > 0. 



Growth of structure 

Models TwCDM la,fb,2a,2b (see Table IJ) were se- 
lected for further study. CosmoMC was rerun with T 
and w fixed, to obtain the best-fit values of the other non- 
derived parameters for input back into CAMB, namely 
fii,/i 2 , ft c h 2 , Hq, n s (scalar spectral index), A s (scalar 
amplitude) and r rc i (optical depth of reionization). 

Figure [3] shows the CMB power spectrum of the best- 
fit parameter sets for the chosen values of T and w. The 
only significant difference between the CMB spectra is in 
the ISW feature, although this is not very large because 
CosmoMC has fit them well to the data from WMAP7. 
By contrast, there are dramatic differences between the 
total matter power spectra at z — for these models 
(see also Fig. 2]). We chose not to fit the matter power 
spectrum to observational data - because the modifica- 
tion to the growth of matter perturbations S m due to 
the interactions is degenerate with the galaxy-DM bias 
b in observations of galaxy number density fluctuations: 
5 g = bS m . This degeneracy is governed by equations (pM)) 
and (|25|). Figure [3] does not include any bias. 

Note that 51 c can be very small in models with large 
r, since it can be compensated for by a higher w in order 
to obtain a sensible Hq. This explains the correlation in 
the r — w plane shown in Fig. [21 so the late-time effect 
of the DM may be proportionately even greater than one 
might think at first glance. 



The combination of similar ISW signatures and large 
differences in the growth of structure is unusual - in a 
ACDM cosmology for example, different growth rates 
lead to correspondingly dissimilar ISW signatures. The 
mechanisms behind this are clearest from the growth 
of DM perturbations in the Newtonian limit: on sub- 
Hubble scales at late-times, 

<5 C »0 = V, S x = (/>' = ^ = 0, (42) 

in the Newtonian gauge (B = = E). The evolution 
of synchronous gauge density perturbations in CAMB 
matches that of perturbations in the Newtonian gauge. 
The ISW effect comes from gravitational potentials de- 
termined by the Poisson equation, 

fc 2 = ~Ai:Ga 2 { Pc 5 c + p b 8 b ), (43) 

and the left panel of Fig. [4] shows that there are indeed 
large differences between the models in their growth rates 
at late times. The reason the ISW effects can remain 
small for these models is that the non-standard back- 
ground evolution (see Fig. |6|) can counteract the growth 
of 5 C in (|4"3"|) and lead to relatively stable gravitational po- 
tentials. The right panel of Fig.fJ]shows that the relevant 
combination, a 2 p c 5 c , can remain comparable for models 
with very different structure formation histories such as 
those considered here. Note how well the TwCDM la 
model mimics the ACDM behaviour of a 2 p c S c , effectively 
leading to the same \ 2 (see Table I). 

This important feature of IDE models has implications 
for any cosmological test which assumes a standard evo- 
lution of the DM energy density during matter domina- 
tion, such as those for detecting deviations from GR. It 
may also be useful for distinguishing between IDE and 
modified gravity models |l04| , which have standard back- 
ground evolutions. 

Using (P), 42]), 4321, pj, |34]), g3l) and the Fried- 
mann equation, % 2 — 8irGa 2 ptot/3, a velocity indepen- 



dent equation of motion for 
(5||u c model: 



S r can be derived for the 



l-^^)s^4nGa 2 {p b 6 b 



2p to t aT p x 



2 - 3w 



aT ( 



)]}■ 



(44) 



3p c H p 

Thus the DM perturbations experience effectively differ- 
ent values of % and G due to the interactions: 



Has 
~W 
Gen 
~G~ 



1 



= 1 



a£(>x_ 

2ptot «r p x 

3p c H p c 



2-3w 



£|1 
H 



Px 
Pc 



(45) 



(46) 



The model by contrast has a non-standard Euler 

equation (|36p , and there remains a term proportional to 
9 X which can not in general be neglected: 



6" + n 



1 - 2^-^)6' c = ^Ga 2 { Pb 5 b + Pc 5 c 



2p t pt «r p x 

3p c H p c 



H p 

2-3w 



aT- 

n 



)]} 



(47) 



Nevertheless, for stable models 6 X remains small enough 
to be negligible and we can define the deviations from 
standard growth due to the interactions via 



H cS 

n 

G e s 
~G 



1 



2 ^£Px 

H p c ' 
2p t ot ar p x 

3p c H p c 



aT 

3W+ H 



(48) 
(49) 



These equations show that the differences in mo- 
mentum transfer lead to a greater modification to the 
growth via H e g for the QHitx model and via G e ff for the 
Q\\u c model, as can be seen in Fig. [5] It is clear that DM 
perturbations in the models with large couplings are al- 
ready beginning to grow e xponentially at the present day 
(compare @,|M[9a[l0|]). In models with Q£ = I> c <, 



as studied in 73, SJ, 84|, there is no interaction source 
term in the synchronous gauge version of (|32t and so the 
DM perturbations are stable. 



V. CONCLUSIONS 

We have studied a model of dark sector interactions 
with an energy transfer proportional to the DE energy 
density, and with momentum transfer vanishing either 
in the DM or the DE rest frame. We performed an 
MCMC analysis and found the best-fit parameters us- 
ing a data compilation that predominantly constrains the 
background evolution. We found model constraints to 
which ACDM is a good fit, although parameter degen- 
eracies do allow for significant interaction rates at the 
present day and even admit the two extreme cases of 
zero DM at early times and zero DM today. 

We analyzed the growth of structure in this model and 
found that the effects of large growth rates on the ISW 
signature in the CMB can be suppressed by the non- 
standard background evolution. We also showed that 
interactions can greatly enhance growth in these models 
via effective Hubble and Newton constants, in varying 
degrees depending on the momentum transfer. 

There appears to be some tension between the back- 
ground evolution and structure formation. The CMB, 
SNe and BAO data slightly favour interactions, while the 
growth rate of DM perturbations likely rules out large in- 
teraction rates. There is a degeneracy with galaxy bias, 
which deserves further investigation. This would allow 
the use of full range of large-scale structure data and 
would significantly improve the constraints on the IDE 
models considered here. 

Interacting models are known t o be deg enerate with 
modified gravity models [7lL 173. f!^5l4T07l | . It is impor- 
tant to break this degeneracy, in order to strengthen cos- 
mological tests of GR - currently devised tests do not in- 
corporate the possibility of a dark sector interaction. The 
key distinguishing features of IDE and modified gravity 
(MG) occur in: (1) the late-time anisotropic stress, i.e. 
4> — ip; (2) the evolution of the background DM density, 
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p c (l + z)~ 3 ; (3) the DM-baryon bias: 



<t> ~ "0 

5b - S c 
Ob — 8 C 



MG 
^0 

= const 
= const 
= 



IDE 

= 

7^ const 
7^ const 
can be nonzero 



These features are the basis for breaking the degeneracy. 
For example, any difference between the metric potentials 
can be tested via peculiar velocities (a probe of 4>), weak 
lensing and ISW (both sensitive to <f> + ip). 
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Appendix A: MCMC Analysis 

Using CosmoMC we explore the full parameter space 
of the TwCDM model, as illustrated in Fig. [7] The code 
was modified to vary the two new parameters T and w 
and further coding was necessary to ensure that models 
with negative DM energy densities were rejected from the 
MCMC analysis. In addition, the scaling solution for the 
background used to find the redshift of the two BAO data 
points was replaced by coding to take into account the 
non-standard background evolution of the models. Cos- 
moMC varies a parameter 9 — 100 times the ratio of the 
sound horizon to the angular diameter distance, in place 
of Hq, because it is more efficient. However the deriva- 
tion of 9 also assumes a standard background evolution. 
We therefore chose not to use 9, but to constrain Hq di- 
rectly instead. Note that the HST prior on Hq assumes a 
particular model of ACDM for evolving H(z = 0.04) up 
to the present day and so has slight model dependence, 
which we neglect. 

There is a plane of degeneracy in the T — w — il c param- 
eter space which allows for an entire range of possibilities 
from zero DM at early times to zero DM at the present 
day - see Figs. [7]and[8l 

The cosmological parameters of the median and best- 
fit models from CosmoMC for w = —1 and when the 



entire parameter space is considered are shown in Ta- 
ble Ull For Q\\u c models, the ISW creates a preference 
in the mean likelihood function for T < 0, as was found 
previously for the Q% = Tp c u^ models [84j . 

For Q\\u x models, there is a preference for T > 0. De- 
spite this the median samples have T < and w < — 1 for 
both the Q\\u c and Q\\u x models. The background data 
therefore shows a slight preference for values of T < and 
w < —1. Both the best-fits and the median samples how- 
ever are relatively close to ACDM, given the wide range 
of interaction strengths allowed. The w = — 1 results are 
included here to show the proximity of V to for these 
models, in line with ACDM. 

The singularity in the perturbations at w = — 1 leads 
us to impose |1 + k;| < 0.01, allowing us to explore the 
entire parameter space. The effect of the w 7^ — 1 insta- 
bility (J3SJ) is illustrated in Fig.[9j The wedged gaps in the 
distribution of accepted MCMC chain steps are given by 
the boundaries of the instability region, defined by (|4"Tj) . 



Appendix B: Initial Conditions 

In synchronous gauge, <f) = B = and ordinarily the 
residual gauge freedom is eliminated by setting 9 C = 0. 
For the Q\\u x model the interaction term in the DM Euler 
equation (f36l) does not in general allow for 9 C = 0. How- 
ever, since T ~ Hq <C H in (|3"2")) - (pT7| . the interactions can 
be neglected at early times. Using 3tp' + k 2 E' = —h/2, 
where h is the synchronous gauge variable [108l | , the evo- 
lution equations used to find the initial conditions for the 
dark sector are, 



2^ + ^ = 0, 6^ = 0, 

5' x + 3H(1 ~ w)5 x + (1 + w)9 x 



(Bl) 



+ 9H 2 (1- W 2 )^ + (1 + W )^ = 0, (B2) 



9' x - 2U9 X - 



k 2 8. T 



(1 + to) 



= 0. 



(B3) 



The dominant growing mode solution for h found in |108| 
leads to the standard adiabatic initial conditions for DM, 

5d - -~h = -\c(kT) 2 ., 6 ci = 0. (B4) 

For DE, we find the leading order solutions, in agreement 
with [Toil, 



C(l + w)k 2 r 2 
12w - 14 



12w - 14 



(B5) 



[1] E. Komatsuet al. (WMAP), Astrophys. J. Suppl. 192, [2] P. J. E. Peebles, AIP Conf. Proc. 1241, 175 (2010), 

18 (2011), 1001.4538. 
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FIG. 7: 68% and 95% contours of the marginalised probability distribution for the Q\\u c model (left) and QHii^ model (right). 
The dashed lines cross at the position of ACDM, the crosses indicate the best-fits in each case and the circles indicate the 
median samples. Note that some areas appear only due to smoothing of the distributions (see Fig. 
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FIG. 8: 68% and 95% contours of the marginalised probability distribution in the Q m — w plane for the Q||u c model (left) and 
Q\\u x model (right). 
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Model 


Qa 


Ax 2 


T/Ho 


w 


Ho 


n b h 2 


n c h 2 


n s 


A s 


Trci 


w = -l TCDM best-fit 


Q\\u c 


-0.146 


0.154 


-1 


70.8 


0.0222 


0.0974 


0.959 


2.16xl0 -9 


0.0824 


w = -l TCDM best-fit 


Q\\u x 


-0.0522 


0.0916 


-1 


70.0 


0.0222 


0.105 


0.959 


2.18xl0" 9 


0.0852 


all T, all w best-fit 


Q\\u c 


-0.294 


-0.806 


-1.23 


70.5 


0.0222 


0.180 


0.956 


2.17xl0~ 9 


0.0822 


all r, all w best-fit 


Q\\u x 


-0.0879 


0.302 


-0.951 


70.1 


0.0222 


0.0823 


0.959 


2.18xl0" 9 


0.0879 


all r, all w median 


Q\\u c 




-1.01 


-1.29 


70.7 


0.0221 


0.194 


0.956 


2.18xl0~ 9 


0.0841 


all r, all w median 


Q\\u x 




-0.578 


-1.19 


70.7 


0.0222 


0.164 


0.958 


2.18xl0" 9 


0.0840 



TABLE II: Cosmological parameters of the median and best-fit samples from CosmoMC for w = — 1 and when the entire 
parameter space is considered. 
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